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Abstract 

The intrinsic decay rate of orthopositronium (o-Ps) formed in Si02 powder was measured using a mod- 
ified method which determined the time dependence of the pick-off annihilation rate using high-energy- 
resolution germanium detectors. That is, the main systematic error due to thermalization was accounted 
for and integrated into the time spectrum fitting procedure. The averaged value was found to be 
A _ Ps = 7.03991° ™^ (stat.) ± 0.0015(sys.)^s~ 1 , which agrees well with the 0(a 2 ) QED prediction; vary- 
ing 2.8-4.2 experimental standard deviations from other measurements. 



1 Introduction 

1.1 o-Ps lifetime problem 

The bound state of an electron and positron, positro- 
nium (Ps), is a pure quantum electrodynamical sys- 
tem that provides a highly sensitive field for test- 
ing accurate descriptions of bound states in quan- 
tum field theory. The triplet (l 3 Si) state of Ps, or- 
thopositronium (o-Ps), will in most cases decay into 
three photons due to having odd-parity under the C- 
transformation. For the last two decades, vast num- 
ber of theoretical works @, § § § §, |, |, § have been 
done to yield the decay rate of o-Ps, up to some par- 
tial results on 0(a 2 ) correction. Recent theoretical 
efforts have evaluated the value up to 0(a 2 ) correc- 
tions (7.039 934(10) ^s" 1 jl), order a 3 ln(l/a) cor- 
rections (7.039 968(10) fis^ fR], O), and so far ne- 
glected binding energy corrections JL2| which is apart 
from P] by about 70 ppm. Therefore there still ex- 



ists a 100 ppm level uncertainty on the theoretical 
side. The decay lifetime is sufficient such that a large 
number of direct measurements have been carried out 
|l3| , [Til |l5| , Jitj , 17 , Il8| . However recent precise mea- 
surements |19[ 20, |2l[] have indicated that the mea- 
sured decay rate is much larger, i.e., 1670, 1600, and 
1170 ppm respectively, than the theoretical QED pre- 
diction |J. These respective discrepancies of 9.1, 8.0, 
and 5.2 experimental standard deviations are statisti- 
cally significant, being referred to as the 'o-Ps lifetime 
problem.' 

Although a variety of experiments have since been 
carried out to search for exotic decay modes of o- 
Ps @ H El HI' EH |7|, Hi, none have provided 
confirming evidence to elucidate this discrepancy. 

1.2 Thermalization of o-Ps 

We went on to perform further study [ p9| in which 
two independent precision measurements were car- 
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ried out using two types of SiC>2 powder, 
correlated well, giving a combined value of 



Results 



A _ Ps = 7.0398 ± 0.0025(stai.) ± 0.0015(sys.) (is' 1 , 

which is consistent with the QED 0(a 2 ) predic- 
tion. This measurement employed a completely new 
method that took into account the thermalization 
process of o-Ps. 

Since some fraction of o-Ps inevitably results in 
'pick-off annihilations due to collisions with atomic 
electrons of the target material, the observed o-Ps 
decay rate A b s is the sum of the intrinsic o-Ps decay 
rate A3 7 and the pick-off annihilation rate X P i C k, i-c, 



X bs(t) — A3 7 + \ p ick(t), 



(1) 



where X P i C k is proportional to the collision rate of 
o-Ps with the target materials; hence being propor- 
tional to the densities of the target materials and the 
velocity of o-Ps. 

Immediately after formation, o-Ps has a kinetic en- 
ergy of about 1 eV, thermalizing via elastic collisions 
with surrounding molecules to an ultimate thermal 
energy of about 0.03 eV. Due to the deceleration or 
thermalization process, X P i C k must be described as a 
function of time dependent on the properties of the 
respective materials. A fc s was previously measured 



9, 20, 21fl by varying the densities of the target ma- 



terials, with the extrapolation to zero densities ex- 
pected to yield A3 7 under the assumption of constant 
o-Ps velocity. Unfortunately, however, such extrapo- 
lation results in large systematic error due to an ex- 
pected low thermalization rate, especially regarding 
the low density limit. Measurements in SiC>2 pow- 
der indicated this systematic problem ||(| ; while 
in other experiments, o-Ps thermalization rates mea- 
sured in gas were also substantially smaller than pre- 
viously thought (3l). Another problem is that the 
assumption of linear dependence of X p i C k on the den- 
sities of the target materials is not always valid, i.e., 
o-Ps collisions with materials could very well be more 
complex due to the occurrence of multiple scattering 
at material surfaces. 

Therefore, our previous measurements adopted a 
formulation containing all information representing 



the thermalization process [£9| ; hence completely re- 
moving systematic errors related to the thermaliza- 
tion process. From Eq. (|l|), the population of o-Ps, 
N(t), can be expressed as 



N(t) = N exp -A 



Xpick (t ) 



A 



3 7 



(2) 

Note that the energy distribution of photons from the 
3-body 7-decay is continuous below the steep edge at 
511 keV, whereas the final state of the pick-off anni- 
hilation is 2-body; a behavior producing the 5ll-keV 
monochromatic peak. If the energy and timing in- 
formation are simultaneously measured, \ p i C k(t) I A37 
can be determined from the energy spectrum of the 
emitted photon. Then, providing the ratio is deter- 
mined as a function of time, the intrinsic o-Ps de- 
cay rate A3 7 can be directly obtained by fitting the 
observed time spectrum with Eq. (^); an approach 
which precludes any ambiguity due to the thermal- 
ization process. 



2 Experimental setup 

2.1 Uncertainties in previous mea- 
surements 

Although our 1995 results were consistent with QED 
predictions, the error was large in comparison with 
that of other experiments ^0[ ^l| . The time spec- 
trum was obtained using CsI(Tl) scintillators which 
provided a relatively broad time response function 
(<r~4 ns) and long scintillation decay time. Accord- 
ingly, since some fraction of the prompt annihilation 
may have flown into the o-Ps decay region, uncer- 
tainties occurred in the time spectrum range imme- 
diately after the prompt events. This uncertainty led 
a restriction in the statistics such that we obtained a 
rather large error. 

Here, although a new experimental setup/tech- 
nique is used similar to the previous one, the above- 
mentioned problems are mitigated by employing 
Nal(Tl) versus CsI(Tl) scintillators. 
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p=0.03g/cm 3 
Figure 1: Experimental setup. 

2.2 Apparatus 

Figure [l] shows the experimental setup, where a very 
weak 22 Na (3 positron source (0.08 fiCi) with 2-mm 
spot diameter is sandwiched between two sheet-type 
plastic scintillators (NE104) having a 12-mm diam- 
eter and 100-fim thickness. The source, held by a 
cone made of 25-/im-thick aluminized myler, is sit- 
uated at the center of a 48-mm diameter vacuum 
chamber made of 500-^m-thick glass. The chamber, 
evacuated to 4 x 10~ 2 torr, is filled with SiC>2 sil- 
ica powder, AEROSIL 300CF (obtained from Nippon 
Aerosil Ltd.), having a grain size of 7-nm in diame- 
ter with density of 0.03 gjcrr? . Since the powder is 
hydrophile, they were heated for 4 hours at 180°C to 
remove the absorbed water molecule. 

Most emitted positrons passing through the scin- 
tillators transmit a light pulse to a trigger photo- 
multiplier (PMT) (Hamamatsu H-3165-04) such that 
Ps is formed as the positrons are stopped in the sil- 
ica powder. The photons subsequently emitted by 
the decay of o-Ps are measured using two types of 
detectors. 

Two germanium semiconductor detectors (diame- 
ter, 61 mm; thickness, 69 mm; Ortec, GEM 38195) 
precisely determine A P i C fc(i)/A3 7 as a function of time. 



The detectors are arranged perpendicular to each 
other (Fig. to reduce the incidence of Compton 
scattered photons from the other detector. A 2.0- 
mm-thick lead sheet located in front of the detectors 
eliminates the contribution of the two simultaneous 
low-energy photons produced by the 37-decay of o- 
Ps. Due to the detector's excellent energy resolution, 
i.e. ,1.78 keV FWHM at 1274.6 keV, the resultant 
energy spectrum allows the monochromatic 27 pick- 
off annihilation to be easily distinguished from the 
continuous 37-decay. 

The four Nal(Tl) scintillators (crystal size, 
64x64x64 mm) simultaneously obtain time and en- 
ergy information. Because they have much higher 
efficiencies compared to the Ge detectors, while also 
having good time response functions, the obtained 
time spectrum can be fit to Eq. (0). 

2.3 Electronics 

Trigger PMT output is fed to a discriminator which 
provides the start signal for both time-to-digital con- 
verters (TDCs) (Hoshin C006 and CAEN C414) and 
the reset signals for other electronics blocked (vetoed) 
during measurements. The TDCs are calibrated by a 
time-calibrator (Ortec 462) having an intrinsic time 
accuracy of 50 ppm. 

One output from the Ge detector preamplifier is 
fed into a fast-filter amplifier (FEA) (Ortec 579) whose 
output is divided into three signals by a linear fan- 
out module (LFF)(Lecroy 428F). One signal feeds a 
constant-fraction discriminator (CFD) (Ortec 473 A) 
whose output is used as the stop signal for the TDC, 
while the other two feed a charge-sensitive analogue- 
to-digital converter (C-ADC) (Lecroy 2249SG). This 
enables measuring energy with a narrow gate width 
(1.25 /is) and the base-line condition just prior to 
an event having a 200-ns ADC gate width. The dif- 
ferential and integral time of the fast-filter amplifier 
are optimized to obtain good time resolution of 6 ns 
rms. To measure a precise energy spectrum, the other 
output from the preamplifier is amplified by a spec- 
troscopy amplifier (Ortec 673) such that shaped sig- 
nals are fed into a peak-holding amplitude-to-digital 
converter (PH-ADC) (Hoshin C011) with 40- /is gate 
width. 
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The output from each of the four Nal(Tl) scintil- 
lators is divided by the LFF into four signals. One 
output is fed into a discriminator whose output with 
good time resolution of 1.2 ns rms provides the stop 
signal for the TDCs, while the other three outputs 
are fed to C-ADCs. One C-ADC with gate 
width, called the wide ADC (Lecroy 2249 W), mea- 
sures the whole charge for the duration of the signal; 
one with 250-ns gate width, called the narrow ADC 
(Lecroy 2249SG), measures the charge itself; and one 
with 180-ns gate width, called the base ADC (Lecroy 
2249SG) measures the base-line condition of the sig- 
nal just prior to the event. 

Accordingly, each Ge detector and Nal(Tl) scintil- 
lator provides three types of ADC information. Be- 
cause pile-up events result in larger charge integration 
on the wide versus narrow ADC, their values are com- 
pared such that these events are effectively rejected. 
The base ADC is also used to reject pile-up events 
which can occur if the base-line condition is shifted 
due to any remaining tail of the previous signal. All 
ADCs use entirely independent gate input such that 
they are able to start charge integration using self 
timing, which results in the ADC gate width being 
narrower in accordance with signal duration. 

Approximately the o-Ps events were collected for 
two month, and during the time a room temperature 
was maintained at 26.0±0.5°C to ensure stability of 
the amplifiers, ADCs, and TDCs. 

The absolute strength of magnetic field around the 
Ps formation assembly was measured to be 0.5 ± 
0.1 Gauss corresponding to a mixing between the o- 
Ps and the p-Ps by the ratio 3 x 10~ n . Thus the 
observed o-Ps can be considered as a pure o-Ps sam- 
ple. 

3 Analysis 

Use of the energy spectrum precisely measured by the 
Ge detectors enables determining \ p i C k(t)/ \ 3j . The 
energy spectrum of the photon originating from the 
pure o-Ps sample is obtained by subtracting the acci- 
dental contribution (time range, 2000 — 3000ns) from 
the measured spectrum. Figure |^ shows the resul- 
tant o-Ps spectrum, where the pick-off annihilation 
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Figure 2: Energy spectrum of o-Ps decay obtained 
by Ge detectors using a time window of 160 — 710 
ns. The solid line indicates the 37-decay spectrum 
calculated by Monte Carlo simulation, agreeing well 
with measured data below 508 keV. 

peak at 511 keV is clearly evident at the edge of the 
37-decay continuous spectrum. 

Using Monte Carlo simulation in which the setup 
geometry and various material distributions are re- 
produced in detail allowed us to calculate the ex- 
pected 37-decay energy spectrum. That is, for every 
simulated event, three photons are generated accord- 
ing to an 0(a) corrected energy spectrum, with suc- 
cessive photoelectric, incoherent, or coherent scatter- 
ing interactions of every photon with materials being 
followed until all photon energy is deposited. The re- 
sponse function of the detector is determined based 
on the measured spectrum of monochromatic 7-rays 
emitted from 152 Eu(344.3 keV), 85 Sr(514.0 keV), and 
137 Cs(661.7 keV), the lines of which are convoluted 
in the simulation. 

As shown in Fig. |[ the o-Ps and 37 spectra indi- 
cate good agreement below 508 keV. Although the 
simulation produced a detailed structure of the actual 
experimental setup, the shape of the 37 spectrum is 
nearly independent of simulation details in that the 
sharp drop-off at 511 keV is solely due to the phase- 
space cutoff, being almost entirely determined by the 
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Figure 3: Enlarged view of measured spectrum after 
subtracting the 37 contribution from the o-Ps spec- 
trum. 



energy response function of the Ge detectors. Figure 
|3| shows an enlarged view of the o-Ps spectrum af- 
ter subtracting the 37 spectrum, where the resultant 
spectrum peak and measured Ge detector response 
function fit well. Since the center of the fitted peak 
is located at 510.86 ± 0.21 keV, this indicates that 
subtraction of the 37 contribution is properly per- 
formed. Hence, the peak is considered to represent a 
pure sample of pick-off annihilations, and the number 
of pick-off annihilations can be counted. The number 
of decays of o-Ps into 37's can accordingly be counted 
using the o-Ps spectrum such that X P i C k is sim- 
ply estimated from the counts ratio, where n^^jn^ 
is multiplied by a normalization factor considering 
relative efficiencies of the Ge detectors and the 0(a) 
corrected energy spectrum. Since \ P i C k/ ^37 is only 
dependent on the relative efficiencies of the Ge de- 
tectors vice absolute efficiencies, the resultant ratio 
is stable against various systematic errors. 

ratios are respectively calculated in the 
same manner using various time windows, with Fig. 
U showing time dependencies of these values. Be- 
cause the fractional energy loss of o-Ps per collision 
is dependent on its energy in Si02 powder, the time 
dependence of the average kinetic energy of o-Ps at 
time t, E(t), can be expressed by the following dif- 



time (ns) 

Figure 4: Pick-off annihilation rate and 37 decay rate 
plotted as a function of time. Fitting results (solid 
line) are obtained using Eq. The horizontal line 
at Xpick/M-y = 0.787 x 10~ 2 shows the approximate 
value at t — » 00. 

ferential equation |}2], ^3| : 

d- 

dt 



E(t) 



2m Ps E(t) 



' K ' 3=0 



(3) 

i/2 



where mp s ,Al, and L are the mass of o-Ps, effec- 
tive mass at the surface of the Si02 grain, and mean 
distance between the grains. Since the pick-off anni- 
hilation rate is proportional to the average velocity 
of o-Ps, Api C fe/A3 7 (= 9(t)) can be expressed as 



dt 



6(t) = -C(6{tf 



E 

i'=0,l,2,- 



a 3 9{t)\ (4) 



where #oo = 6{t — > 00), C, and a,j(J — 1,2, • • •) are 
constants. With the exception that the summation 
term is replaced with abbreviated form 9(t) K in which 
k is an arbitrary real number, Eq. (||) is used for 
fitting measured Xpick/X^ values. As indicated by 
the results in Fig. |^, o-Ps takes about 600 ns to 
become well thermalized in Si02 powder. 
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Figure 5: Time spectrum of Nal(Tl) scintillators with 
a photon energy window of 350-450 keV. 

Figure |5| shows the time spectrum of Nal(Tl) scin- 
tillators using a photon energy window of 350 — 450 
keV, where the sharp peak from prompt annihila- 
tions is followed by exponential decay of o-Ps and 
the flat spectrum of accidental events. Due to using 
a weak positron source (0.08 (J>Ci), the o-Ps decay 
curve covers a wide region of about 1.2 ^is , being 
about eight times as long as the o-Ps lifetime. 

Using the least square method the time spectrum 
is fitted with 



F(t) = exp(-R stop t) 



1 



tpick ^pick 
£37 ^3 7 



N(t) + C 



(5) 

where N(t) is in Eq. (||), and e V i C k and £3 7 are re- 
spectively the probabilities that the photons emitted 
from the pick-off annihilation and 37-decay of o-Ps 
deposit energy of 350 — 450 keV in the Nal(Tl) scintil- 
lators. Results of Monte Carlo simulation estimated 
£picfe/e37 to be 0.128 ± 0.002, while R s t op is the mea- 
sured stop rate of the Nal(Tl) scintillators. 

Open squares in Fig. || show fitted values of the 
decay rate of o-Ps, A 37 , plotted as a function of fitting 
start time using a simple exponential function, i.e., 

/(*) = exp(-i? stop ) [ No exp (-A o6s t) + C]. (6) 
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Figure 6: Decay rate of o-Ps plotted as a function 
of fitting start time. Open squares and closed circles 
show values obtained by fitting with Eq. (j^) and the 
proposed method, respectively, while the horizontal 
line shows the value obtained at the start of fitting 
(140 ns). 



Since observed A Q fc s contains the pick-off rate fraction, 
the intrinsic decay rate A3 7 is obtained by correcting 
Xobs by a factor of (l + X P i C k(t — > oo)/A3 7 ) « 1.00787. 
Closed circles in Fig. || show fitting results obtained 
using the proposed method, where the obtained val- 
ues are stable against the fitting start time. Figure 
[j] shows the effect due to thermalization has been 
completely taken into account in our new method. 
The horizontal error-bars indicate statistical errors 
solely due to fitting the Nal(Tl) time spectrum, while 
the vertically extending bars represent propagated er- 
rors due to the fitting of Api C fc(i)/A 37 . Due to the 
tail effect of prompt events, fitting x 2 is rapidly in- 
creasing before 120 ns, instead one can obtain sta- 
ble x 2 values after 140 ns. Note the stable fitting 
values after 140 ns, where the fitting \ 2 1S 621.308 
for (573-3) degrees of freedom; hence this time is 
used as the fitting start time. The resultant value is 
A 37 = 7.0401 



■-0.0033 
-0.0038 



(stat.) [is , which includes the 
statistical error in determining X p i C k{t) / X^ . 

When systematic errors were estimated for vari- 



G 



cms effects, the predominate contribution was found 
to be due the Monte Carlo simulation. That is, 
the 37 contribution in the o-Ps spectrum is esti- 
mated by normalizing the Monte Carlo distribution 
with measured data from 480 — 505 keV. In this 
measurement, we used the same Monte Carlo code 
adopted previously 29 . Equivalent systematic er- 



= 7moi+r^i(stat.) + _ o o ° \l- 



rors originating from the simulation are accordingly 
expected and estimated to be +190 and —150 ppm. 
Other errors arise due to in-uniformity of Si02 pow- 
der (±30 ppm), fluctuations in relative efficiencies of 
Ge detectors dependent on the position of o-Ps for- 
mation (—40 ppm), and the position of o-Ps decay 
(—50 ppm). The remaining sources of error are due 
to hardware considerations: the absolute accuracy of 
the time calibration (±50 ppm), uncertainties in the 
relative efficiencies of the Ge detectors (±50 ppm), 
dependence on the base-line condition of Nal scintil- 
lators (±23 ppm), and relative efficiency of the Nal 
scintillators (e P i C fe/e3 7 ) (±7 ppm). Since these sys- 
tematic errors are independent of each other, the to- 
tal systematic error is their quadratic sum. 
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Figure 7: History of the measurement of o-Ps decay 
rate. The solid line shows 0(a 2 ) corrected QED pre- 
diction, while the dashed line shows 0(a) corrected 
value. 



resun is a 37 - _ 003B 

. Total error is obtained by the quadratic 



The result is 
(sys.) fis 

sum of them, i.e., A 37 = 7.0401+°;™^. Although 
the obtained precision is relatively lower compared 
to previous experiments, the results are essentially 
consistent as shown in Fig. [7] which indicates the 
history of o-Ps decay rate measurements. In addi- 
tion, substantial improvement in the timing system 
is demonstrated. 



4 Conclusions 

The o-Ps decay rate was measured in Si02 powder 
using a modified method in which the effect of ther- 
malization of o-Ps is accounted for by high energy- 
resolution Ge detectors. By combining obtained and 
our previous result p9[ that assume systematic er- 
rors mainly stem from the same origins, we obtained 



A 



37 



7.0399+-- 



° (stat.) ± 0.0015(s?/s.)^s- 1 . This 
is consistent with the 0(a 2 ) corrected QED predic- 
tion of 7.039934(10)Ais" 1 [§ within an error of 360 
ppm, while varying 2.7 ~ 4.2 a from other recent 
measurements [[l9| |2(j, |2lJ . Realizing that this exper- 
iment represents dynamic development, further im- 
provement in Ge timing response is expected soon; 
an innovation that should enable us to investigate 
the fast component preceeding a 140-ns start time 
such that the 0(a 2 ) contribution in QED can be ul- 
timately be confirmed at a level of just 200 ppm. 

Sincere gratitude is extended to Professor Toshio 
Hyodo and Dr. Yasuyuki Nagashima, University of 
Tokyo, for their insightful suggestions and valuable 
discussions concerning o-Ps thcrmalization. 
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